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CL(' Abstract. Axford and McKenzie [1992] suggested that the energy released 

in impulsive reconnection events generates high frequency Alfven waves. 

The kinetic equation for spectral energy density of waves is derived in 

the random phase approximation. Solving this equation we find the wave 

J>-^, spectrum with the power law "-1" in the low frequency range which is 

i-G [ matched to the spectrum above the spectral brake with the power low "- 

1.6." The heating rate of solar wind protons due to the dissipation of Alfven 

waves is obtained. 



1. Introduction 
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0^ I The solar wind which is formed by the gas-dynamic expansion of solar 

^^ ■ corona into the interplanetary space can be divided into two states: the fast 

^ , (> 700 km/s) and slow (~ 400 km/s) streams. The fast solar wind originates 

in the polar coronal holes, regions of the solar corona with relatively low 

^.^. electron temperature T ~ 10^ K, reduced plasma density n = 10® cm~^ 

^S^ \ and magnetic field at the base of corona -Bq ~ 10 G [1], [2]. Variations of the 

Mh' velocity of this fast solar wind streams are insignificant (700-800 km/s). The 

slow solar wind is associated with transient openings of closed field regions 

in the corona. The streams of slow wind are limited by the range ±13° 

5^ ■ near the equator. The transition zone between this regions has the width 

5^ I ±13° — ±20° on the north and south latitudes in the solar minimum [3]. 

The radiation balance analysis in the near equatorial coronal holes [4] shows 

that to accelerate solar wind it is necessary to dissipate the energy flux of 

radiationless origin « (5 ± 1) • 10^ erg/cm^ -s within the distance of one- 
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two solar radii above the Sun's surface in order to maintain the observed 
temperature 1,5- 10^ K of slowly rising expanding gas. 

Axford and McKenzie [1] suggested that the energy source, necessary 
to accelerate fast solar wind, is in the regions of strong magnetic field 
that define the boundaries of chromospheric supergranules. If the magnetic 
field is not strictly unipolar in this regions, for example, the closed loops of 
magnetic field lines are existing at the base of the corona, then the necessary 
energy is released in the processes of impulsive reconnection of magnetic 
field lines on the characteristic spatial scale of network activity / = 100 km. 
Such processes of impulsive reconnection are accompanied by the generation 
of Alfven and fast magnetosonic waves. Using the parameters of plasma at 
the coronal base described above we can find the Alfven wave velocity Va = 
B / y/A^Tmrnl = 2 • 10^ cm/s, the ion's thermal velocity Vti = 2 • 10^ cm/s, 
and the characteristic time of magnetic field lines reconnection [5] : 

^Va ^ 7r-2-10S cm/s 
^^ 4/ In Rem 4-107 cm • ln3.5 • 10" "" ^' ^' 

where: Rcm = 47rne^VA//w-e^eiC^ = 3.5 • lO^"*^ is the magnetic Reinold's 
number, v^i = 20n/(T°)^' ^ = 2.3 s~^ is the collision frequency. It is evident 
that the characteristic period of Alfven waves is of the order of obtained 
reconnection time. The total energy flux necessary to accelerate fast solar 
wind is estimated as 8 • 10^ erg/cm^ -s [1] and as in the case of the near 
equatorial coronal hole the major part of the energy of generated waves 
should be dissipated within one-two solar radii above the solar surface. 

As the dominant mechanism of dissipation of Alfven waves generated 
by the processes of impulsive reconnection of magnetic field lines, we will 
consider the induced scattering of Alfven waves by plasma ions [6]. For 
the sake of simplicity as in the paper [6] we will limit ourselves by the 
case of the circularly polarized Alfven waves propagating along the radial 
magnetic field Bq = BQ{r)z from the polar coronal holes with the wave 
vector k = (0,0, fc), where r is the radial distance from the center of the 
Sun. With the help of the above cited parameters of plasma and magnetic 
field at the base of the corona, we fix the following ordering of the ion 
cyclotron LOd, Alfven a;^ and ion Doppler kvTi frequencies: 

ujci :^ uJk^ kVri, (2) 

where ujd = eBo{r)/miC, ujk = kVA- Under this conditions plasma is magne- 
tized and the resonant condition describing the interaction of Alfven waves 
with the different polarization propagating in the opposite directions along 
the magnetic field lines can be written as: 

^k + ^k' - {k + k') Vz = 0, (3) 
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where u>k = ^Va > 0, cuk' = —^'Va < 0, k,k' > and v^ is the velocity along 
the magnetic field. Here in the process of induced scattering the quantum 
of higher frequency is exchanged for that of lower frequency. As a result of 
this some fraction of wave energy is transferred to plasma ions. 

In the next section we will derive the kinetic equation for the interact- 
ing waves in the random phase approximation. In this derivation we limit 
ourselves by the case Tj » Te [7], [8], that is valid for the solar corona. 
Then in the section 3 we will describe the evolution of the Alfven waves 
spectrum in the processes of their propagation and show that the wave en- 
ergy flux necessary to accelerate the fast solar wind can be dissipated at 
the heights of 1-2 solar radii above the solar surface that is required by the 
Withbroe [4] analysis. 

2. Derivation of kinetic equation in the random phase approxi- 
mation 

The kinetic equation for the particle distribution function fj in the absence 
of collisions takes the form [9], [10]: 



d d d \ „ iei ^4. ^jg 



^ + ^-TJT ^ ^-J7^ • fi = ^^.^k(^^ 



dt dz '86 J ■" 2mjC 

dvz \ k ^ ) v±d9 



Wfc ^.\ ^ , 



k J dv± 

■ fj exp [-i {ujkt - kz)] 



(4) 



where: v±_ and 9 are the velocity and the asimuthal angle of the cyclotron 
rotation of particles, B^ and B^ are the waves amplitudes with the left 
and the right polarizations propagating away from Sun and towards the 
Sun respectively. 

Following perturbation theory we solve the kinetic equation (4) by find- 
ing the distribution function of plasma particles as power series in the wave 
amplitudes: 

fj = fjo + fji + fj2 + fj3, (5) 

where fjQ is the isotropic maxwellian distribution of particles of j species 
in the zero approximation and /,„ is the particle distribution function on 
the n-th step of iteration. Using the iteration procedure for solving the 
kinetic equation (4) we obtain the corresponding Fourier's components of 
the particle distribution functions for the first two steps of iteration: 
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w, 



/,2 = -^Va 



+ 



{k + k')VA {k^ + k'^)Vl jk-k') 



^CJ 



2a;2, ' 2uo. 



UJ 



cj 



dfjo/dv. 



C] 



ujk + ujj./ — [k + k') Vz + iO 



(8) 



^2^' exp [-i (lja; + Wfc') t + i (/c + A;') z] . 



Here we have neglected the sniah contribution related to the differentiation 
of nonresonant denominator on Vz in the equation (8). 

Perturbation of the charge density related to the particle distribution 
function of the second order fj2 excites the longitudinal wave with the 
potential 4>k+k'i that can be found from the Poisson equation: 



{k + fc')^ (pk+k'^k+k' = 4vr ^ Sj / (fvfj2, 



(9) 



where ek+k' is the dielectric permeability of the plasma: 



efc+fc' 



1 + E 



w, 



p] 



{k + k') 



d^v 



{k + k') dfjo/dvz 



LOk + ujk' — {k + k') Vz + iO 



(10) 



With the help of the operator of particle interaction with the longitudinal 
waves in a form: 



/. 



{k + k') (pk+k' dfjo i(u^^+u;,/)t+i(k+k')z (-.-.^ 



rrij ujk + ook' — {k + k') Vz + iO dv. 



we find on the next iteration step the particle distribution function of the 
third order for the waves propagating from Sun: 



/; 



J3 



2m jC 



dk' 



^k' 
k' 



d _d_ 

dv_L dvz 



2tt 



uJk - kvz 



iOcj + ^0 



{fj2 + fj4 



(12) 



xB,*e *^ exp [i{ujk't — k'z)] . 



The kinetic equations for the waves propagating from Sun and towards 
the Sun respectively can be described in terms of Maxwell equations for 
the Fourier's components of the wave amplitudes: 



Bt 



^^EeW^v.,e-^(/J. + /3- 



(13) 



Here we take into account that the second order current, associated with 
the distribution functions fj2 and fj^, does not contribute to (13). 
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In the linear approximation the dispersion relation (13) can be written 
as [11]: 

Taking the frequencies ordering (2) we obtain from the equation (14) the 
well-known dispersion relation for the Alfven waves: 

In the order to obtain the time dependent evolution of Alfven waves we 
multiply the equation (14) by B^*{k,cd)eyip[i{oj — to)] and integrate it by 
dCbduj. Using the relation e^^^ = e^^ + ie^^ , we obtain: 



jj duduel^B+{k,u:)B^ 



~^ i(^-^)t 



+ V,,.\^|R+|2 (16) 



"2 du^k dt +'^ki^^k)\B^\^ 

where the time dependent wave amplitude is defined as: 

BUt) = J du;B+{k,co)e-^^--^-^\ (17) 

Restoring the contributions of the nonlinear currents on the right-hand side 
of the equation (13) we find the kinetic equation for the weakly interacting 
waves propagating from the Sun: 



dt ^-^ kc^ J 271 J ^k — kvz — ujcj + iO ng") 

To obtain the above equation in such form we integrate by parts the expres- 
sions for nonlinear currents to get rid of the derivatives over Vz and v± and 
neglect small corrections related to the differentiation of the nonresonant 
denominator by Vz- 

Expanding the expressions for the nonresonant denominator and for the 
particle distribution functions fj2 and fj^f^ as a power series in cok/i^cj and 
LOy /iOcj^ kvz and k'vz to the second order we find the equation for waves 
propagating from the Sun in the form: 
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dt ' '= ' ^ 2c2 i 27r y k \ k' J \ ujh uj. 
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cj 



dfj'Q/dv 
x{ l + A 



E ^l^: I d^V (1 + Av) — - 






r 



uJk + uJk' - {k + k') Vz + iO 



df,o/dvz \B^\'\BZ\' 



UJk + UJk' — {k + k') Vz + iO Bq 
where 



(19) 



_ {k + k')VA , {k^ + k'^)Vi , {k-k')Vz , {k'^-k^)VAVz 
^j ~ o, , . + o, ,2 + o, , . + , ,2 



2wci 2u;^j. 2u;ci ^cj 

and Aji and differ from Aj only by the index f . 

Considering the hmit Tg <C Tj we obtain in the last term in figure 
brackets that the electron contribution to integrals dominates that for ions: 



E/^^v 



dfj'o/dvz 



UJk + ^k' — (k + k') Vz + iO 

^pe f ,3 9fjo ^' 



' r Bf uj^ ^'^^ 

pe I ,3 'JjjO ,, '^pi 

' a V — :^ — » 



k + k' J Vzdv-z (k + k') vj^^ ' 



As the result the first term in the brackets cancels with the last one ap- 
proximately equal —1. Taking the integral by parts over v^ and then over 
dk' we reduce equation (19) to the equation obtained by Livshits and 
Tsytovich [7], [8]: 



J^\Bk I - -'^'=^2-^1^. I k. (21) 



In the case of stationary expansion of plasma from the solar corona along 
approximately radial magnetic field lines this equation takes the form: 

Hr) + V'.Ml llBtf = -kvJ^^JBilH. (22) 

3. Evolution of Alfven wave spectrum in the processes of their 
propagation 

Evolution of Alfven wave amplitudes in the processes of their propagation 
into the interplanetary space we calculate theoretically assuming that (a) 
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the interplanetary magnetic field BqIv) and the solar wind velocity u are 
radial and spherically symmetric, (b) the Alfven waves are circularly po- 
larized, (c) wave amplitudes are small \B^{r)\'^k <C BQ^r). 

Integrating the equation for the spectral energy density of waves^ 



[««^-.(")i:"''r' 


fc ...2y.jBUrrkd\B 


we find the general solution [12]: 


\B+{r)\'k = ^ 


, j VA{r)dr 1 1 
J u + VAir) k{RQ) k 

Rq 



{r)\^k 



Att 



(23) 



(24) 



Here: $ is an arbitrary function; VA{r) = VA{RQ)RQ/r; k{RQ)VA{RQ) /^ir = 
1.7 Hz; Rq = 6.96- 10^° cm; n = |S^(r)p/|B+(r)|2 is the ratio of the spec- 
tral energy density of Alfven waves, propagating to the Sun and from the 
Sun respectively, that is the constant because both waves have the same 
spectral index a = —1 and radial profiles. The ratio of the wave energy 
density to the energy of the magnetic field / = \B^{r)\'^k/BQ{r) happens 
to be the constant too due to the identical radial profiles of the energy 
density of Alfven waves and the energy of magnetic field and very weak 
dependence of the wave energy density |i?^(r)p/c on wave vector k due to 
the very small left-hand side of the equation (23). Numerical value for / 
can be found by equating the energy flux of Alfven waves generated at the 
base of the corona to the energy flux of the radiationless nature with the 
value 8 • 10^ erg/cm^ -s, necessary for the acceleration of the fast solar wind 
from polar coronal holes: 

IB+Pyfc 
VA(r) ' ^' =8-10^ erg/cm^-s. (25) 

47r 

As the result we find that: / = 5 • 10~^. 

In the process of wave propagation to the interplanetary space the po- 
sition of the spectral break, separating spectral domains with the different 
slopes, is shifted to the low frequencies [13]. From the general solution (24) 
we construct the particular solution: 

J__.f yA{r)dr _J_ 

hr~ J u + VA{r)^ kiR^y ^^^^ 

Rq 

^In the works [7], [8], authors were interested in stationary spectra of Alfven waves 
when the left-hand side of the equation has been zero, in the other words d/dt = 0. 
That's why they didn't call any attention that the sign in the right-hand side of the 
equation should be mines not plus. 
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We can find the exact position of tlie spectral break frequency as a function 
of radial distance by dividing both sides of this equation by the Alfven wave 
velocity: 



fbrir) 



0.6 + 2ttkI 



Rr. 



VA{r)dr 
u + VA{r) 



Rp, 



(27) 



Assuming that the acceleration of fast solar wind is terminated at r = 
3ii!0 [16] we find that fbri'^Re,) ~ 0.28 Hz. At this stage most of the heating 
has been deposited within (1 — 2)Rq above the surface as required by 
Withbroe [4] analysis. We can simplify this expression at the distance r > 
WRq assuming that the velocity of the fast solar wind u = 7.5 • 10^ cm/s 
^ VA{r). As a result we reduce the expression (27) to the form: 



fbrir) 



0.6 + ttkIu 



Rp) 



R, 



■0 



(28) 



Taking the numerical value of the spectral break 6 • 10~^ Hz at the distance 
r = 6Oi?0 [13], separating the domain of Alfven wave turbulence with the 
spectral indexes a = —1 and domain with a = —1.6 and equating it to 
the theoretical value obtained from the equation (28) we find the ratio of 
the spectral energy density of Alfven waves propagating to the Sun and 
from the Sun respectively that is equal to k = 0.17. This means that the 
intensity of waves propagating to the Sun is six times lower than those for 
the waves propagating from the Sun. 

Experimental observations show that above the spectral break a wave 
spectrum with a slope near —1.6 [14] is established: 



\Bt^ 



r)p 



\B2. 



{r)\'klf/k'-'. 



(29) 



The equation for the spectral energy density of waves in this case should 
inherit the structure of the equation (23) with the continuos transition 
across the spectral break: 



[u{r) + VAir)] 



dr Att 

-Kk'-'^VA{r)k^; 



o.2\B^{rrk'-' d \Bt{r)?k'-' (^0) 



Bl{r) 



dk 



4tt 



Equation for characteristics of (30) take a form: 



1 dk _ ^j[27TfUrT^VA{r)''-^ 



A;i-6 dr 



u+VA{r) 



(31) 
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Figure 1. The characteristics of the equation (30). 



where / = \B^^^{r)\'^hr/ Bl{r) = 5 • 10-^ u > VA{r). Beyond the 187.5i?0 
(0.87 astronomical units — AU) the collisional hydromagnetic turbulence 
spectrum of Kraichnan [15] with the spectral index a = —1.5 [13] replace 
the collisionless spectra with the spectral indexes a = —1 and a = —1.6. 

Integrating (31) and transforming the result to the frequency depen- 
dence on the radial distance we obtain: 



C 



Re 



■& 



0.6 



1 



fbr{r) 



0.6 



(32) 



where C is an arbitrary constant. 

Figure 1 shows the characteristics of the equation (30) between the 
curves of the spectral break frequency f^r and the ion cyclotron frequency 
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fci = 1.6 • 10*^ {Rq/v) , where the wave energy is totally absorbed in the 
ion cyclotron resonance. 

To estimate the fraction of Alfven wave energy transferred to the protons 
of solar wind plasma in the process of induced scattering of waves we use, 
following Tu [17], the equation for the magnetic moment of protons. Doing 
this we take into account that due to the smallness of the left-hand side 
of the equation (30) the product \B'^{r)\'^k^'^ has a very weak dependence 
on the wave vector k. Therefore the right-hand side of (30) also depends 
weakly on k and we can rewrite it in form: 

Hr) + VA{r)] |:^^ = -^^Q.2f>KfUr)l'Bl{r) =. -^n(r), (33) 

where n(r) is the volumetric energy flux. 

The equation for the magnetic moment take a form: 

d ^ fT^\ n(r) 

where Tx_ is the perpendicular proton temperature in Kelvin, ks Boltz- 
mann's constant, n{r) = IQ^^Rq/t)'^ cm~^ density of plasma. 

We integrate the equation (34) from the three solar radii assuming that 
the acceleration of solar wind is already achieved. As a result we have: 



^- ^- ' "(•■' .* = ai441n 



Bq Bq J n{r)kBuBo{r) 

3R0 



J. 
0.6 + ttkIu . 



(35) 



Non-adiabatic heating of the protons related to the presence of the 
volumetric energy flux n(r) reach only 0.16 K/nT in the range of radial 
distances [3-^0, 200/?©]. Therefore the source of energy for the acceleration 
of fast solar wind from the polar coronal holes must now be found in the 
solar corona itself [18]. 

Below the spectral brake frequency fbrif) we have the spectrum of 
Livshits-Tsytovich [7], [8] with the spectral index a = —1. However this 
spectrum doesn't extend from /f,^ to the whole low frequency domain but 
truncated by the WKB solution. 

Solving the stationary equation for Alfven waves [19], [20], [21]: 

V ((Snir) ± 2V^(r)) ^^^^) - u(r)V^?^^ = 0, (36) 

the equation for the wave amplitude was found in the form: 



u 



n(r)-3/2|5+(r)|2A; = const. (37) 
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At the distances r > 10i?Q we can neglect the small parameter VA{r)/u <C 1 
and rewrite it in the form defining the reduction of the wave amplitudes 
and their spectral index a = —2/3 [13]: 

const • n{r) = (\B+{r)\^kf'^ . (38) 

Besides that, knowing from observations that the spectral break curve fbr 
and the WKB solution are merging at the distance 0.87 AU (187.5 Rq) 
we are able to find the form and position of WKB solution curve from the 
equation (38) in the form (fig. 1): 

A™M = 1.6. 10-3. (^^)'". (3,,) 



4. Conclusion 

Using the coUisionless kinetic equation for interacting Alfven waves in the 
random phase approximation we have followed analytically the evolution of 
these waves in the process of their propagation from the polar solar corona. 
As a result of such evolution the boundary of waves with the spectral index 
a = — 1 is shifted towards low frequencies due to induced scattering of 
waves by ions, which increase energy of ions by the recoil effect. However, 
this spectrum does not extend over the whole low frequency range but is 
truncated by WKB solution. 

Above the frequency fbr the spectrum with the spectral index a = —1.6 
is obtained. Let us note that both in the polar solar corona and in the 
interplanetary space the approximation of coUisionless plasma is valid at 
least up to the distance 0.87 AU. 

A. A. Galeev thanks the Alexander von Humboldt Stiftung for support 
and Prof. Joachim Triimper for attention to this work. This paper was par- 
tially supported by the RFBR grants 97-02-16489 and 96-15-96723 (science 
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